Noise-assisted quantum transport and computation 
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The transmission of an excitation along a spin chain can be hindered by the presence of small 
fixed imperfections that create trapping regions where the excitation may get caught (Anderson 
localization). A certain degree of noise, ensuing from the interaction with a thermal bath, allows to 
overcome localization (noise-assisted transport). In this paper we investigate the relation between 
noise-assisted transport and (quantum) computation. In particular we prove that noise does as- 
sist classical computation on a quantum computing device but hinders the possibility of creating 
entanglement. 



I. INTRODUCTION 

Decoherence and relaxation induced on a quantum sys- 
tem by a fluctuating environment can considerably affect 
its dynamical behavior. The quantum transport of a par- 
ticle in a one dimensional tight-binding model, for exam- 
ple, requires coherent hoppings between nearest neigh- 
bor locations. The presence of imperfections appearing 
as a random time-independent potential can hinder the 
transmission. In fact, the energy mismatches due to such 
imperfections lead to destructive interference of the wave 
ftmction fT] and Anderson localization appears. 
However, it has been recently argued that, in certain sit- 
uations, the presence of an environment can assist the 
transport of an excitation The authors consider 

exciton transport on small complexes, such as the Fenna- 
Matthew-Olson (FMO) protein complex and engineered 
quantum systems (binary trees) showing that, at room 
temperature, a limited amount of decoherence (account- 
ing for the interaction between the complexes and the 
phonons of the surrounding environment) cancels out the 
localizing effect of small random imperfections along the 
transmission line. See also @. 

In this note, we apply similar considerations to the model 
of universal quantum computer put forward by Feynman 
in one of the earliest paper on quantum computation Q . 
This computer model suits very well the investigation of 
the relation between transport and computation, since 
the transmission of an excitation does correspond to the 
execution of a computation. 

The device consists of an array of two level systems 
with nearest-neighbor interactions mediated by ancillary 
qubits. The ancillae play the role of an input/output reg- 
ister. An excitation traveling along the chain plays the 
role of a cursor administering the applications of unitary 
transformations (computational primitives) to the regis- 
ter. In particular, when the cursor reaches the end of the 



chain, the register is in the output state and the compu- 
tation is complete. 

We investigate the effects of irregularities in the on-site 
energies of the cursor chain and the possibility of assisting 
the transmission of the cursor along the chain by means 
of an external potential and of energy exchanges with the 
environment. We show that the cursor does indeed travel 
beyond the localization length and reaches the far end 
of the chain. We discuss the conditions on the system- 
environment interaction and on the kind of information 
processed that allows for an extension of noise-assisted 
transport to noise-assisted computation. 

The paper is organized as follows: we first discuss the 
detrimental effect of random static noise on transmis- 
sion (section ini and present an instance of noise-assisted 
transport by means of energy exchanges (section lllB . We 
then discuss the effects of noise on the computational 
capabilities of the Feynamn quantum computer (section 
IIV[) . The last section is reserved to discussion and out- 
look. 



II. LOCALIZATION ON A LINEAR LATTICE 

Consider an excitation moving along a one dimensional 
regular lattice, e.g. an electron moving along a periodic 
potential, in the tight-binding approximation. Calling 
I X ) the state localized at site x, the Hamiltonian for a 
lattice of s sites is: 



1 

x=l 



The evolution | "0* ) of an initial condition | i/jo ) is: 
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where: 
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, fc = 1,2, ...,s 



sm I — — I I .T , 



s + 1 



(2) 
(3) 



We are interested in the transport capabihty of the chain, 
that is in the probabihty that an excitation initiahy lo- 
cated at the left hand of the chain (| V'o )=| 1 )) reaches 
the rightmost site s, namely in the probability 

PQ = s{t)^\{s\ijt)\\ 

where Q is the position operator 



distribution with mean zero and variance a^. With this 
choice, the localization length is [l^ 

The idea of adding a linear potential Vl = 
~9Y^x=i 2^1 2; )( a; I, with g > 0, to reduce the effect of 
the irregularities in the chain comes from a classical prej- 
udice: classically, one expects this potential to pull the 
excitation rightward; moreover, with a suitable choice of 
the magnitude of g, the defects (small random potential 
terms) would become negligible with respect to the linear 
potential. 

It is immediate to see that this idea/prejudice is very 
naive and actually quite misleading: the eigenstates of 
the Hamiltonian 



x=l 

Figure [T] shows an example of the motion of the exci- 
tation along the lattice: the particle moves ballistically 
along the chain until it bounces against the boundaries. 
As time goes on, the wave packet spreads over the line. 
The probability of finding the excitation at the end of the 
chain is highest when the first reflection occurs (around 
time t — s) and, at that time, is of order 0{s~'^^^) [1]. 

We now turn our attention to the effects of imperfections 
on transmission. We model imperfections of the lattice 
by a random time-independent potential. In the position 
representation we adopt, the random potential appears 
in the main diagonal of the Hamiltonian: 



x){x\ 



(4) 



The coefficients e^, x — l,2,...,s are independent real- 
izations of a random variable R. 

The fact that the presence of imperfections can local- 
ize the excitation/cursor (Anderson localization) is well 
known [ol-fl^. The phenomenon can be understood in 
different ways. One is in terms of destruction of the co- 
herent hopping between nearest neighbors, that leads to 
destructive interference of the wave function. Another 
interpretation, that better suits our discussion, is given 
in terms of the eigenstates of the Hamiltonian Hn: the 
presence of small imperfections on the chain results in a 
localization of the eigenstates of Hr- In our particular 
instance, since the initial condition has a non-negligible 
projection on the eigenstates localized at the left of the 
chain, and the eigenstates do not change shape (but only 
phase) in time, the excitation remains localized around 
its starting position. 

The typical size of the localized states (localization 
length) depends on the distribution of the noise. Here 
we consider the random variable R to have a Gaussian 



Hr+l =Ho + Vr + Vl = 

s 

= Ho + ^{ex - gx)\x){x\ 



(5) 



with g > are still localized. Localization manifests 
itself, in this case, in the form of Bloch oscillations, a 
genuinely quantum effect: the momentum of the quasi- 
particle moving on the lattice changes linearly with time, 
p{t) = p{0) — gt, until it reaches the boundary of the Bril- 
louin zone, where it is Bragg reflected [l^. The localiza- 
tion len gth is determined by the intensity of the force 
through [in : 
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where Ah^. is the bandwidth of the Hamiltonian Hr. 
As an example, figure [5] shows distribution of the observ- 
able Q for the eigenstates of Hr+l, for s = 20, ct = 0.5 
and 5 = 2: the eigenstates are localized. If the particle is 
initially localized at site 1 it has a non-vanishing projec- 
tion only on the eigenstates localized in the same region, 
and will therefore remain localized there. 
The eigenvalues of the same Hamiltonian for a sample of 
100 realization of the random imperfections are shown in 
figure db). 

In the presence of imperfection, therefore, the transmis- 
sion of a particle along a linear lattice via Hamiltonian 
evolution is suppressed and the idea of pulling the par- 
ticle rightward by means of an external constant force, 
mimicking the effect of a battery, does not work. 



III. NOISE ASSISTED TRANSPORT 

The onset of interaction with a thermal environment 
might be expected to further reduce the transmission ca- 
pabilities, of the system. But this is not necessarily the 
case. Recent work showed that some interaction with the 
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environment can assist transport on short chains [3| and 
quantum networks 0, 0] ■ 

In this section we investigate the effects of dissipation on 
the transport capabilities of the linear lattice described 
in section HIl 

The linear chain is now weakly coupled to a thermal bath. 
The full Hamiltonian (system+environment) is: 



H' 



TOT 



Hr 



(6) 



where Hs is the Hamiltonian of the chain (either Hq or 
Hfi or Hfi^ L ) , He is the Hamiltonian of the environment 
and Hi describes the system/environment interaction. 
Following [l^, we suppose Hj = —£ ® S, where £ is 
an observable of the environment and S is an observable 
of the system. 

In order to know the specific form of £ and S we should 
know the precise implementation of the array and the 
nature of the environmental degrees of freedom. Since 
this is not the case, we limit ourselves to formulating the 
hypotheses under which we derive the results of this and 
the following section. 

• The environment must act as a reservoir: the ef- 
fects on it of the interaction with the system are 
negligible. Moreover, we assume that the reservoir 
is in a thermal stationary state, that is oc e~^^'^ 
at some inverse temperature /3. 

• The system exchanges energy with the environ- 
ment, not particles. In other words, the number 
of excitation in the chain remains a constant of the 
motion even when it interacts with the reservoir. If 
the operator S is of the form: 



where the dissipator V has the form: 



m— 1 n—1 

— —L{m, n)^ L{m, n)p — —pL{m, n) f L{m, n). 



Here C is the (positive) system-environment cou- 
pling constant. The generators L(m, n) are defined 
in terms of the normalized eigenstates of Hs (satis- 
fying: Hs\n) = e„|n); ei < e2 < ... < e^; {n\m) = 

L(n,n)=0 (10) 
L{m, n) = I m ) ( n I, for m ^ n. 

Each coefficient ^{m,n) in ([S]), for m ^ has 
the meaning of transition probability per unit time 
from \ n) to | to ). We set, in what follows, for 
n ^ ra, Lu(m,n) = u>{n, m) = |e„ — em|- 
We require, following Il7|, that transitions n ^ m, 
with n > TO, in which the system absorbs energy 
from the reservoir, take place at the absorption 
rate: 



7(to, n) — t] {uj{m, n)) ■ 



1 



,{l3 ui{m,n)) _ ^ ' 



(11) 



Transitions n ^ m, with to < n, in which the 
reservoir drains energy from the system, take place 
at a rate 



5 = ^/(:r)|x)(:r:| 



(7) 



x = l 



for example, this requirement is satisfied. We point 
out that this choice is not unique but is, by far, the 
simplest, since it assumes that the bath interacts 
locally with the chain sites. In the following section 
this requirement will have significant implications 
on the computational capabilities of the system as 
well. 

• The environment can possess a continuous spec- 
trum, but due to the weak coupling the system 
couples only with the resonant levels. Moreover 
we require the density of states of the bath to be 
constant over the resonant energy levels. 

• The evolution of a given initial state p(0) of the 
system is determined by the Lindblad equation: 



dpjt) 
dt 



^t[p{t),Hs]+V{p{t)), 



(8) 



7(m, n) ^ Tj [uj{m, n)) 



1 , (12) 



taking into account both stimulated emission and 
spontaneous emission, via the additional term +1 
appearing in (fT2|) . 

• We consider only the case of an underdamped mo- 
tion in which simultaneous emission of two or more 
phonons is forbidden; namely we assume that 

r]{uj{m, n)) = d,n.n+l + Sm,n-1, 

but for irrelevant constants that we absorb in the 
system/environment coupling constant C. 

Under these conditions, the equations for the elements of 
the system density matrix p{t), in the energy representa- 
tion, are the following: 
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-T.P^iAi) = -i{e{m) - e{n)) - 



j{j,m) +"t{j,n) 



dt 
d 

di' 



p,n,n{t), m^n 



(13) 
(14) 



c=l 



where e(k), k 
tonian Hs [H 



1 , . . . , s are the eigenvalues of the Hamil- 



If the Hamihonian of the system is the free Hamiltonian 
Hq, the transmission capabihties of the computing de- 
vice is dramatically reduced by the presence of the bath, 
as expected: the system relaxes toward a thermal state 
that, the lower the temperature 1//3, the more is concen- 
trated around the middle of the chain. 
The introduction of disorder is obviously not expected to 
do any better: the eigenstates of the Hamiltonian are 
localized around the peaks of the random potential. 
The situation is radically different when a linear (ex- 
ternally tunable) potential is added to the random one: 
Bloch localization makes the eigenstate of the Hamilto- 
nian Hji+L localized in space (see figure [2ja)). The in- 
teraction with the bath determines a "stroll" over the 
energy eigenstates which, because of the structure of the 
bath, favors, at low temperatures, the propagation to- 
ward states of lower energy. The localized nature of the 
eigenstates of i^i^+L, due, we remind, to the presence of 
the external potential Vj,, makes the random walk on the 
energies quasi-equivalent to a random walk on the sites of 
the chain: losing energy amounts to moving to the right. 

The speed of the walk on the energy landscape induced 
by the dissipation depends on the energy gaps between 
the eigenstates of -fffl+L and on the coupling constant (. 
This dependence can be analytically accessed by inspec- 
tion of equations [13] and [131 

The diagonal elements describe a random walk on the 
energy landscape of the system. The presence of a suf- 
ficiently strong external potential {g > 1) makes the en- 
ergy gaps |e(fc) — e{k ± 1)| almost equal for fc = 1, . . . , s 
(see, for example figure[2Ub)); the emission (and absorp- 
tion) rates 'jsix) = j{x, x + 1) (and 7^(2;) = 7(x -I- 1, x)) 
are then independent of x (7_B(a;) « jE^lAix) ~ 7a)- 
For /? high enough, becomes negligible and 7^; « 1. 
Therefore the random walk on energies described by 
equation [TH proceeds toward states of lower energy at 
a rate C 7-E ~ C- 

The coherences (off-diagonal terms) are described by au- 
tonomous equations (|T^ and are exponentially depressed 
in time (decoherence). 

Frames (a) and (b) of figure [3] show the evolution of the 
initial condition po = I 1 )( 1 I under the Lindbladian 



r 



(|5]) for P — 1, g = 2 and a realization of the random 
imperfections extracted from a zero-mean Gaussian with 
standard deviation a — 0.5. 

The transport process determined by the Hamiltonian 
dynamics (see figure [TJ is qualitatively different from the 
one depicted in figure [3] 

Under "perfect" Hamiltonian evolution the particle 
bounces back and forth. It is thus necessary to deter- 
mine the optimal time at which to measure the position 
of the excitation, that is the time at which the probability 
of finding the excitation at the end of the chain is max- 
imal and this probability decreases with s as 0(s~^/^). 
Moreover, the presence of disorder suppresses transmis- 
sion. 

In the presence of a bath and an external static field, the 
transmission of the particle to the far end of the chain 
takes a time that is one order of magnitude larger than 
under Hamiltonian evolution in the absence of imperfec- 
tions (compare the time scales figures [1] and [3]) . This is 
due to the weak system-bath coupling = 0.05). But 
once the excitation has arrived at the end of the chain it 
remains there, and the motion is stable with respect to 
random fiuctuations of the on-site energies. 
The coupling with a cold reservoir, in the presence of a 
uniform force field that completely offsets a weak ran- 
dom potential, assists the transport (however slow and 
incoherent) of an excitation along a line. 



IV. NOISE-ASSISTED (QUANTUM) 
COMPUTATION 



In this section we discuss the relation between quan- 
tum transport and quantum computation. The quantum 
computer model we use in our investigation is the one 
proposed by Feynman Q . 

Feynman's quantum computer consists of two logically 
separated parts; one part, the clock, is an excitation mov- 
ing along a lattice. The second part, the input /output 
register, is a collection of additional degrees of freedom, 
say n spin 1/2 particles cr(j) = {(Ji{j),cr2{j),cr3{j)), 3 = 
l,2,...,n. 

The system is governed by the time-independent Hamil- 
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tonian: 



1 '"^ 

Hp ^--'^\x + l){x\^U^ + \x){x + l\ ^U;^ 



(15) 



The presence of a (random) potential acting on the cur- 
sor can drastically reduce the probabiUty of finding the 
excitation at the end of the chain (i.e. the computation 
completed) but does not alter the working mechanism of 
the Feynman machine. 

For the sake of definiteness, let the state | ipo ) evolve 
under the Hamiltonian 



Each term of the Hamiltonian involves two nearest neigh- 
bor sites of the clock and a unitary operator Ux acting on 
the register. The ordered product Ug-i . . .U2U1 realizes 
some input/output transformation we want the quantum 
computing device to accomplish. Figure S] shows the ar- 
chitecture of the machine. 

Because of the properties of the Hamiltonian Hp, the 
position of the excitation along the chain uniquely deter- 
mines the state of the register. This fact has interesting 
consequences. Consider an initial condition of the form 
I -00 ) = I 1 ) I i?i ) = I 1, i?(l) ), i.e. with the particle 
located at the beginning of the chain and the register in 
an input state | i?i ). Then the set 



H 



F,V 



V = Hp 



s(0o) = {|i,i?i),i2,i?2: 



S,Rs)}, 



where \ Rj ) = Uj-i . . .U2Ui\ Ri ), constitutes an or- 
thonormal basis, the computational basis, or Peres basis 
, for the region of Hilbert space visited by the evolved 
state \ ipt) = exp{—iHpt)\ ipo ). We refer to the space 
spanned by the Peres basis as to the computational sub- 
space. 

In particular, if upon measurement the cursor is found 
at the rightmost site of the chain, the register collapses 
to the output state | R{s) ) = Us-^i . . . U2Ui\ R{1) ). 

It is clear that the capability of the chain of transfer- 
ring the excitation from one end to the other is a central 
matter in the analysis of the computational power of the 
Feynman machine. 

Before discussing the kinematics of the cursor, we point 
out that the sole effect of the interaction of the clock with 
the register (of n spins) is the appearance of a degeneracy 
of order 2" in the spectrum {ek}l^i of the tight-binding 
(clock) Hamiltonian Hq (see ([T])). Once an initial con- 
dition of the form | i/jq ) = I Ij -R(l) ) has been set, how- 
ever, the spectrum of the Hamiltonian Hp restricted to 
the computational subspace is no longer degenerate. In 
this subspace, to each eigenvalue ek (see ^) there cor- 
responds the eigenvector 



\Vk) = 



1 ^ 

x=l 



sm 



klTX 



X, R{x) 



(16) 



x=l 



f{x)\ x){x 



(17) 



where V represents a potential. The projector on the 
space spanned by the Peres basis B{'ipo) 

s 

P(H^J = ^|x,i?,)(a;,i?, I 

x=l 

remains a constant of the motion. Indeed 



(18) 



X 

= J2fi^)i\^){^l\y,Rv){y,Ry\] = 

x,y 

^J2f{x){\x){x\y,Ry){y,Ry\ 
-\y,Ry){y,Ry\x){x\=0). 



The study of the evolution of the Feynman machine can 
therefore be reduced to the analysis carried out in section 
ITTl In particular, the Hamiltonian ([T7| . restricted to the 
subspace "H^q, and represented in the same basis, admits 
s eigenvalues ek and corresponding eigenvectors | Cfe ). 

By requiring the operators Ux to act on a single regis- 
ter qubit, the Hamiltonian ([15]) becomes 3-local. This 
restriction simplifies the implementation of the device; 
but at the same time it reduces the computational ca- 
pability of the Feynman machine: it does not allow the 
realization of any two-qubit gates. In order to restore the 
universality of the model Feynman introduced the Switch 
circuit [7], that implements the selection statement. For 
a recent review, see p^ . 

In what follows, we consider a particular instance of the 
Switch: the Controlled NOT (CNOT) gate. The gate 
acts on two register qubits: cr(c), the controlling qubit 
and cr(p), the controlled (or passive) qubit. In our im- 
plementation, if the controlling qubit is "up" , namely 
CT3{c) = -1-1, the passive qubit a{p) is "negated" by apply- 
ing CTi(p); if the controlling qubit is "down" (0-3(0) = —1) 
then the controlled qubit is left unchanged. 
The Hamiltonian realizing the gate is: 



6 



He NOT {a)^-]^[\a + l){a\® i±|^ + |a + 2)(a+l|®ai(p) + |a + 5)(a + 2|® i±|^^ (19) 
Ma + 3)(a|® ^^ + |a + 4)(a + 3| + |a + 5)(a + 4|(8) j + H.c. 



Here a indicates the site of the clock at which the switch 
starts and H.c. stands for Hermitian conjugate. 
The state of the controhing qubit determines which 
branch of the CNOT circuit is visited: if a{c) is "up", 
the excitation will walk along the upper path and (Ji(j>) 
is applied. If ct(c) is "down", the excitation will walk 
along the lower path and no operator is applied to the 
register. 

In order to discuss the effects of imperfections and dissi- 
pation on such a circuit, it is expedient of put some "iner- 
tial" sites before and after the CNOT, "during" which we 
do not apply any transformation to the register qubits. 
The full Hamiltonian of the circuit is: 

a— 1 
x=l 

+ HcNOTia) 
1 

"2 E + + + 

x—a+b 

Figure [5] shows the structure of the CNOT gate. 
If the initial state of the system is 

\i'o) = \Q = ^)(^\ ^3(c) - +i,'j3{p) = -1 > = 

= +1,-1). 

only the upper branch of the CNOT is visited during 
the evolution of the system and, if the cursor reaches the 
region a + 5, . . . , s the state of the register is | +1, +1 ). 
The analysis of the motion of the cursor is once more 
simplified by the use of the Peres basis 

B{ij^)={\l)®\ +1,-1),..., 

\a + l)(g)\ +l,-l),|a + 2)(8)| +1,+1), 
|6)® I +l,+l),...,|s)(»| +1,+1)). 

We indicated by 6 = a + 5 the position at which the 
CNOT ends. 

A similar argument applies to the initial condition 
|^^) = |1)(^| -1,-1). With this initial condition 
it is the lower branch of the CNOT that is visited. The 
corresponding computational subspace is spanned by 

e(V>o^)=(|l)0|l,_l),..., 

|a + 3)® I -l,-l),|a + 4)®| -1,-1), 
|6)(8| -l,-l),...,|s)(8| -1,-1)). 



Both bases consist of s — 2 orthogonal elements. The 
Hamiltonian ([20| restricted to H^u (or T~L^d) and rep- 
resented in the basis B {^iPq){oi B {ip^)) is equivalent 
to the "free" Hamiltonian |l]) on s — 2 sites. In par- 
ticular, the motion of the cursor along the two com- 
putational paths is exactly the same. The presence of 
random imperfections, namely of independent realization 
£1, e2j • ■ • J £s of a zero- mean Gaussian random variable R 
with variance ct^, perturbs the motion along the com- 
putational paths. Since ea+i 7^ £0+3 and ea+2 7^ £0+4 
with probability 1, the spectra of the Hamiltonians re- 
stricted to 'H^u and "H^d will be slightly different. But 
the correctness of the computation is guaranteed by the 
conservation law (ITSl) . 

The problem is the ensuing localization that we discussed 
in section II; it can suppress the probability of reaching 
the region 6, . . . , s. 

In the previous section we showed that we can make the 
excitation travel beyond the localization length, and ac- 
tually reach the rightmost sites of a linear lattice, by 
means of an external force and dissipation. In what fol- 
lows we show that this result extends straightaway to 
each computational path. 

We add an external field such that the potential differ- 
ence between one site and its right next-neighbor (s) is 
-ff, by setting 

/a+2 s \ 

Vi^-\y^9x\^){^\+ E 9{x~2)\x){x\\. 

\a:— 1 x— a+3 / 

The Hamiltonian evolution of the system is given by 

s 

HcNOT + Vr + VI^ Hcnot + ^ 6x1 a; ) ( a: I (21) 

x=l 

Ca-2 s \ 

J2gx\x){x\+ J2 9{x~2)\x){x\\. 
x—l x—a-]-3 / 

The initial condition of the system (either | ip^ ) or | ipQ ) ) 
selects the relevant computational subspace (either H^c/ 
or T-L^d). By switching on an interaction with a thermal 
bath satisfying the hypotheses listed in section IIII[ we 
determine a walk on the energy landscape of the selected 
space identical to the one we described in the previous 
section. For example, if the initial state of the machine 
is Pq' = I V'o^ ) ( V'o' 1 1 ttis elements of the density matrix 
(t), restricted to the computational subspace H^u and 
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there expressed in the energy representation, are given by 

dni) and dm). 

This means that the clocking particle can travel beyond 
the CNOT. And when it does, the conservation law (fT8|) 
guarantees the state of the register is the correct one. 
We can thus state that noise can assist classical compu- 
tation on a quantum computing device. 
Figure [B] shows the probability of reaching, starting from 
site 1, the region 6, . . . , s under different evolutions. The 
presence of imperfections along the chain strongly re- 
duces the probability of ever reaching the sites located 
after the gate. For the ballistic evolution there are times 
at which this probability is close to unity. When the evo- 
lution is dissipative, the time required to pass the CNOT 
is one order of magnitude larger than the first "arrival 
time" under Hamiltonian evolution. But the probability 
of finding the computation completed is monotonically 
increasing to one. 

The equivalence between successful transport and suc- 
cessful computation does not extend to the case that 
sees the controlling qubit a{c) in a superposition of "up" 
(ct3(c) = +1) and "down" (cr3(c) = -1). 
Suppose, for instance, that the controlling spin cr(c) 
starts from the state 

I .X ^1 ^ k3(c) = +l) + |g3(c) =-1) 

I c^i(c) = +1 ) = , 

and therefore the complete system starts from the initial 
condition 

I ,U+D . I V'O^ ) + I ^0^ ) 
\% } ^ ■ 

If the evolution were purely Hamiltonian, with identical 
on-site energies in the two branches of figure [Sj appli- 
cation of the CNOT primitive would bring the register 
{a{c),a{p)) to the maximally entangled Bell state [2Qj] 

The presence of random imperfections on the circuit we 
are considering produces two effects: 

• the introduction of random relative phases between 
the two computational paths. Because of the rel- 
ative phases, when the particle is on the right of 
the CNOT the state of the register can be slightly 
different from the target Bell state. But this effect 
can in principle be "bounded away" for small values 
of the variance cr^ ; 

• Anderson localization of the cursor; this is the key 
issue that we discuss below. 

The transmission of the excitation from one end to the 
other of the circuit can be assisted by dissipation. But 



this time dissipation destroys the computation. The den- 
sity matrix of the system can be schematically repre- 
sented as a block matrix: 



PDoit) 


PDu{t) 


^ PUD{t) 


Puu{t) 



(22) 



The blocks pDD{t) and puuit) represent the evolution of 
the projections of the initial state respectively on H^d 
and H^u . The blocks PDu{t) and punit) correspond 
to the coherences between the computational subspaces. 
Lindblad evolution, according to p3)) . destroys all the 
coherences. This does not affect the computation within 
each computational subspace: the decoherence simply 
transforms a quantum walk on the energy landscape into 
a random walk. But the damping of coherences between 
the computational subspaces destroys the desired entan- 
glement; it transforms the state of the register into the 
maximally mixed state 

2 

This is evident from the fact that the von Neumann en- 
tropy of the register for a machine evolving from | ^q^^ ) 
tends to log(2) as f — cx) (solid line in figure[7]). For com- 
parison, the dashed line of figure [7] shows that, starting 
from I tpQ ), the entropy of the register tends to the value 

competing to a pure state. 

For both initial conditions the entropy reaches a maxi- 
mum about the time i when E{Q) is close to s/2; at this 
time yJvar{Q) is large (see figure [Sljb)). 
For the initial condition | t/jq^ ) , the state of the register 
at time t is a mixture of two states: | -I- 1, — 1 ) ( -f 1, —1 | 
(cursor on the upper branch and at the left of the NOT) 
and I +1,-1-1 I (cursor on the upper branch 

and at the right of the NOT). Hence the value log(2) for 
the absolute maximum at time t of the entropy for the 
dashed line of figure [T] 

If the controlling qubit a{c) is in a superposition of "up" 
and "down" , at time i the register is in a mixture of 
three states: | -I- 1,-1 )( -1-1,-1 | (cursor on the upper 
branch and at the left of the NOT) with weight 1/4, 

1 -I- 1, -1-1 )( -1-1, -1-1 I (cursor on the upper branch and at 
the right of the NOT) with weight 1/4 and, this time, 
I — 1,— 1)(— 1,— 1| (cursor on the lower branch) with 
weight 1/2. Hence the value (3/2) log(2) of the maximum 
reached by the solid line. 

For the initial condition | tpQ^^ ) the steep growth of the 
entropy from to log(2) that leads to the plateau around 
time i = 40 is a witness of the disruption of coherences 
between computational subspaces, as it is easy to check 
by direct inspection of the decay rate of the explicit so- 
lution of EqnlT3]for the coherences. 
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V. CONCLUSIONS AND OUTLOOK 

We presented an example of noise- assisted transport 
on a lattice and discussed a set of hypotheses that allow 
for the extension of noise-assisted transport to noise as- 
sisted computation. We showed that it is indeed possible 
to exploit the interaction of the (Feynman) computing 
device with a reservoir to assist classical computation on 
a quantum device. On the other side, the suppression of 



coherences due to the presence of a bath destroys entan- 
glement. 

We are presently exploring the possibility, offered by 
Time-Dependent Density Functional Theory, of simulat- 
ing the position probability distribution of an open quan- 
tum systems with the one of a system undergoing unitary 
propagation under a time dependent potential [21|. The 
interest of this proposal is in the fact (made possible by 
the extension of the point of view of TDDFT to discrete 
systems [H, [2^ and to open discrete systems that 
the same external time-dependent potential may be made 
to act along two branches of the computation, leading to 
completion without loss of phase information. 
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FIG. 1. s = 20; (a) The initial condition \tpo) ~ \ 1) evolving under the Hamiltonian Ho- (b) The expectation E{Q) = {xl)t \Q\^t) 
as a function of time (solid line) together with E{Q) ± ^^var(Q) (dashed lines). 




FIG. 2. s = 20; g = 2; . (a) The probability distributions of the observable Q in some eigenstates for two realizations of 
Hr+l with a = 0.5 and one with a — (no imperfections), (b) The eigenvalues of Hr+l for a sample of 100 realization of the 
random imperfections for a — 0.5. 
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FIG. 3. s = 20; a — 0.5. (a) The initial condition po = | 1 )( 1 | evolving under Hr+l, with g — 2 and (weakly) interacting 
= 0.05) with a reservoir at inverse temperature 13 = 1. (b) The expectation E(Q) = Tr {Qp{t)) as a function of time (solid 
line) together with E{Q) ± ^var{Q) (dashed lines). 
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FIG. 4. The Feynman machine: the s-sites cursor interacts with the register as described in (|15|) . 




FIG. 5. The CNOT described by equation [20l The state of the controlling qubit cr(c) determines which branch of the circuit 
is visited by the moving particle. This determines the transformation applied to the controlled qubit (j{p): <Ji{p) = NOT(p) 
along the upper branch, the identity along the lower branch. 
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FIG. 6. s = 22: with this choice, each computational path consists of s — 2 = 20 sites, thus allowing for comparison with figures 
[T]and[3l a = 9. The probability of reaching, starting from site 1, the region b, . . . ,s under three different evolutions. Dotted 
line: Hamiltonian evolution determined by Hf; dashed line: Hamiltonian evolution under Hf + Vr, for a realization of the 
random potential Vr with a — 0.5. Solid line: the dissipative evolution determined by the Hamiltonian Hf + V[ + Vr for the 
same realization of Vr, g = 2, coupling C, — 0.05 and inverse temperature /3 = 1. 
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FIG. 7. s = 22; a = 9; C = 0.05; /3 = 1. The von Neumann entropy of the register for (t) (dashed fine) and p^^°{t) 
(solid line). The von Neumann entropy for completely mixed of 2- and 3-state systems are represented as horizontal lines for 
comparison purposes. 



